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Abstract 

We use rigorous nonequilibrium thermodynamic arguments to establish that (i) the nonequilibrium entropy S{Tq) of any 
system is bounded below by the experimentally (calorimetrically) determined entropy Sexpt(Tb), (ii) S'expt(To) is bounded below 
by the equilibrium or stationary state (such as the supercooled liquid) entropy S i scl(T'o) and consequently (iii) S(To) cannot 
drop below Sscl(To). It then follows that the residual entropy Sr, is bounded below by the extrapolated 5 OX pt(0) > <Sscl(0) at 
absolute zero. These results are very general and applicable to all nonequilibrium systems regardless of how far they are from 
their stationary states. 



I. INTRODUCTION 



Nonequilibrium states like glasses from supercooled liq- 
uids (SCLs) are abundant in Nature, whose entropy S can 
only be estimated by calorimetrically measured entropy 
Sexpt, which can then be extrapolated to absolute zero. 
The extrapolated value Sr at absolute zero is commonly 
known as the residual entropy and is normally found to 
satisfy Sr > 0. In practice, one considers the isobaric 
entropy S(Tq) of the system as a function of the tempera- 
ture To of the surrounding medium; see Fig. [T] The exis- 
tence of Sr was first theoretically demonstrated by Paul- 
ing and Tolman [lj; see also Tolman 0|. In addition, the 
existence of the residual entropy has been demonstrated 
rigorously for a very general spin model by Chow and 
Wu J3|. The residual entropy for glycerol was observed 
by Gibson and Giauque [4j and for ice by Giauque and 
Ashley Q. Pauling [6[ provided the first numerical esti- 
mate for the residual entropy for ice, which was later im- 
proved by Nagle Q- Nagle's numerical estimate has been 
recently verified by simulation [8|, |9( . The numerical sim- 
ulation carried out by Bowles and Speedy [1Q] for glassy 
dimers also supports the existence of a residual entropy. 
For a brief review of the history of the residual entropy, 
see flll - Fl8j | - Thus, it appears that the support in favor of 
the residual entropy, see the curve Glassl in see Fig. [21 is 
quite strong. Its existence also does not violate Nernst's 
postulate, as the latter is applicable only to true equilib- 
rium states with a non-degenerate ground state |19l421| . 
Indeed, many exactly solved statistical mechanical mod- 
els show a non-zero entropy at absolute zero. However, 
as of yet, no experiment can be performed at absolute 
zero to experimentally determine the residual entropy; 
in all cases, some sort of extrapolation is required. This 
point should not be forgotten in the following whenever 
we speak of the residual entropy. Despite the above men- 
tioned support for the reality of the residual entropy, it 
has become a highly debated issue in the literature |22| - 
[28j as discussed by these authors. The reason for the 
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FIG. 1: An isolated system So consisting of the system E 
in a surrounding medium E. The medium and the system 
are characterized by their fields To,Po,... and T(t),P(i), ..., 
respectively, which are different when the two are out of equi- 
librium. 



debate is that the relationship among S(Tq), S cxp t(To) 
and the entropy Sscl(To) of the corresponding station- 
ary state is not well understood, and understanding this 
relationship is the main theme of this work. 

In the following, we will speak of "the equilibrium" 
state associated with a nonequilibrium state as the sta- 
tionary (time-independent) state. Depending on the con- 
text, the equilibrium state may represent a true equilib- 
rium state such as a crystal or a stationary metastable 
state such as the supercooled liquid. A nonequilibrium 
state in this work will always be taken as time-dependent. 
Accordingly, 5 (To) above should be correctly expressed 
as S(T), t), and in some cases can be expressed as a func- 
tion S(T(t)) of the instantaneous temperature T(t) [29 - 
l3l| of the system; see Fig. [TJ In this work, we will not 
be concerned with T(t). Thus, we will simply use S(Tq) 
for the nonequilibrium state, knowing well that the state 
continues to change with time. 

For the purpose of clarity, we will consider supercooled 
liquids and associated nonequilibrium states (glasses) in 
the following, but the arguments are applicable to all 
nonequilibrium states. The supercooled liquid undergoes 
a glass transition over a transition range, see Fig. [2] 
over which the entropy falls rapidly with lowering tem- 
perature Tq. The transition region is controlled by the 



rate of cooling so that the glass is a nonequilibrium state 
[3lL l32j ] . As the irreversibility due to the glass transition 
does not allow for an exact evaluation of the entropy, it 
has been suggested [23, [28| that the entropy decreases by 
an amount almost equal to S'r within the glass transi- 
tion region so that the glass (see Glass2 in Fig. [2j whose 
entropy lies below the supercooled liquid) would have a 
vanishing entropy at absolute zero. It has been shown 
by Goldstein [l_l| that Glass2 results in a violation of the 
second law. It should be stressed that if there is ever any 
conflict between the second law [33J] and any other law 
in physics such as the zeroth or the third law, it is the 
second law that is believed to hold in all cases. One can 
also argue that to confine the glass into a unique basin 
in the energy landscape requires microscopic information 
[I5L |32| . l34j ; hence, the particular glass cannot be consid- 
ered in a macrostate. Oppenheim [351 ] has also raised 
somewhat of a similar objection. 

We have drawn the two entropy curves ( Glass 1 or 
Glass2) in Fig. [5]that emerge out of the entropy curve for 
the equilibrated supercooled liquid for a given r bs in such 
a way that Glassl has its entropy above (so that «Sr > 0) 
and Glass2 below (so that Sr = 0) that of the super- 
cooled liquid. The entropy of Glassl (Glass2) approaches 
that of the equilibrated supercooled liquid entropy from 
above (below) during isothermal (fixed temperature of 
the medium) relaxation; see the two downward vertical 
arrows for Glassl. It is the approach to equilibrium that 
distinguishes the two glasses, Glassl and Glass2. Almost 
all experimental investigations leave open the possibility 
that Glass2 may materialize if the irreversibility is too 
large. Our work clarifies the situation. 

It is abundantly clear from the above discussion that 
there is a need to look at the relationship between var- 
ious entropies in Fig. [2] As is customary, we treat the 
supercooled liquid as an equilibrium state, even though 
it not a true equilibrium state; see above. We proceed 
by following the strict second law inequality diS > 
[29 . [30l [3a . [37| , see Eq. ((4]) , and use it to prove the fol- 
lowing results applicable to all nonequilibrium systems, 
regardless of how close or far they are from their equilib- 
rium state: 




FIG. 2: Schematic behavior of the entropy of the equilibrated, 
i.e. stationary supercooled liquid (solid curve) and two possi- 
ble glasses (Glassl-dotted curve, Glass2-dashed curve) during 
vitrification. The transition region between To g and Tog has 
been exaggerated to highlight the point that the glass tran- 
sition is not a sharp point. For all temperatures To < T) g , 
any nonequilibrium state undergoes isothermal structural re- 
laxation in time towards the supercooled liquid. The entropy 
of the supercooled liquid is shown to extrapolate to zero per 
our assumption, but that of Glassl to a non-zero value and 
of Glass2 to zero at absolute zero. 



The Eq. ([T]) is consistent with Glassl but not with 
Glass2. All experiments on or exact/approximate com- 
putations for nonequilibrium systems must obey the 
strict inequalities in Eqs. ([2][T]) without exception. This 
is the meaning behind the usage of "... rigorous ..." in 
the title. The actual values of the entropy are not rel- 
evant for the aim of this work, which is to find the re- 
lationship among different entropies under vitrification. 
Because of the possibility that the systems may be far 
away from equilibrium such as in a fast quench, where 
the irreversible contributions may not be neglected, our 
results go beyond the previous calorimetric evidence 
[l2l . Ha . fl8l [23j | . The systems we are interested in include 
glasses and imperfect crystals as special cases. However, 
to be specific, we will only consider glasses below. 



1. Various entropies obey the following strict inequal- 
ities 

S(T ) > S* cxp t(To) > S SC l(To) forT <T 0g , (1) 

so that the entropy variation in time has a unique 
direction as shown by the downward arrows in Fig. 
H Thus, S(T ) cannot drop below S S cl(T ) (such 
as Glass2 in Fig. [5]) without violating of the second 

law m. 

2. The experimentally observed non-zero entropy at 
absolute zero in a vitrification process is a strict 
lower bound of the residual entropy of any system: 



S R = 5(0) > 5 expt (0) > Sscl(0). 



(2) 



II. ENTROPY BOUNDS DURING VITRIFICA- 
TION 

The vitrification process we consider is carried out at 
some cooling rate as follows. The temperature of the 
medium is isobarically changed by some small but fixed 
ATq < from the current value to the new value, and 
we wait for (not necessarily fixed) time r b s at the new 
temperature to make an instantaneous measurement on 
the system before changing the temperature again. At 
some temperature Tg g , see Fig. [2] the relaxation time 
Trciax, which continuously increases as the temperature 
is lowered, becomes equal to T b s - Just below To g , the 
structures are not yet frozen; they "freeze" at a lower 



temperature Tog (not too far from To g ) to form an amor- 
phous solid with a viscosity close to 10 13 poise. This 
solid is identified as a glass. The location of both tem- 
peratures depends on the rate of cooling, i.e. on r b s - 
Over the glass transition region between Tog and T) g in 
Fig. [2j the system gradually turns from an equilibrium 
supercooled liquid at or above Xb g into a glass at or be- 
low Tog [Sll [32|, |38| • We overlook the possibility of the 
supercooled liquid ending in a spinodal [39]. It is com- 
monly believed that Sscl(O) will vanish at absolute zero 
(Sscl(O) = 0), as shown in the figure. However, it should 
be emphasized that the actual value of Sscl(0) has no 
relevance for the theorems below. 

We will only consider isobaric cooling (we will not ex- 
plicitly exhibit the pressure in this section), which is the 
most important situation for glasses. The process is car- 
ried out along some path from an initial state A at tem- 
perature T) in the supercooled liquid state which is still 
higher than To g to the state Ao at absolute zero. The 
state Ao depends on the path A— >Ao, which is implicit 
in the following. The change dS between two neighbor- 
ing points along suc h a p ath is dS = d e S + d[S in modern 
notation [H, Hfl HI HE E3, S3 • The component 



d c S(t) = d c Q(t)/T = CpdTo/To 



(3) 



represents the reversible entropy exchange with the 
medium in terms of the exchange heat d c Q(t) (in keep- 
ing with the modern notation) added to the glass by the 
medium at time t to the medium at T) and the heat 
capacity Cp. It also represents the calorimetrically de- 
termined change in the entropy in any process. The com- 
ponent 



d;S>0 



(4) 



represents the irreversible entropy generation within the 
system in the irreversible process, and contains, in addi- 
tion to the contribution from the irreversible heat trans- 
fer with the medium, contributions from all sorts of vis- 
cous dissipation going on within the system and nor mally 
require the use of internal variables [29|, [30, [30, SO S3 ■ 
The equality in Eq. Q holds for a reversible process, 
which we will no longer consider unless stated otherwise. 
A discontinuous change in the entropy is ruled out from 
the continuity of the Gibbs free energy G and the en- 
thalpy H in vitrification proved elsewhere [29(. Thus, 
we only consider a continuous change in the entropy as 
shown by the two glass curves in Fig. [5J 

Theorem 1 The experimentally observed (extrapolated) 
non-zero entropy at absolute zero in a vitrification process 
is a strict lower bound of the residual entropy of any 
system: 

S R = S(0) > S expt (0). 

Proof. We have along A— >Aq 



S(0) = S(T ) + JdJS + JdiS, 

A A 



(5) 



where we have assumed that there is no latent heat in 
the vitrification process. Since the second integral is al- 
ways positive, and since the residual entropy Sr is, by 
definition, the entropy S(0) at absolute zero, we obtain 
the important result 

S R = S(0) > S expt (0) = S(T ) + JC P dT Q /T Q . (6) 

To 

This proves Theorem [T] The integral represents the 
calorimetric contribution. ■ 

The strict forward inequality above clearly establishes 
that the residual entropy at absolute zero must be strictly 
larger than S oxp t(0) in any nonequilibrium process. 

Theorem 2 The calorimetrically measured (extrapo- 
lated) entropy during processes that occur when t {, s < 
T r eiax{To) for any T) < T) s is larger than the supercooled 
liquid entropy at absolutely zero 

S expt (0) > S SCL (0). 

Proof. Let Q c (t) = d e Q(t)/dt be the rate of net heat loss 
by the system. For each temperature interval dT < 
below T) g , we have 



\d o Q\ = Cp\dT \= f 




Q, 



dt < \dQ\ (To) 



r re ia x (T ) 
I 



Q 



dt, 



T <T 



On 



where \dQ\ (T ) > denotes the net heat loss by 
the system to come to equilibrium, i.e. become su- 
percooled liquid during cooling at Tq. For T) > T) g , 
d c Q = d c Q C q(T ) = C PtC(i dT Q . Thus, 

JCpdTo/To > /Cp, cq dTo/To. 

To To 



We thus conclude that 



5expt(0) > S S CL(0). 



(7) 



This proves Theorem [2] ■ 

The strict inequalities above are the result of glass be- 
ing a nonequilibrium state. We have now verified the 
second statement in the Introduction. 

The difference Sr— S cxp t(0) would be larger, more ir- 
reversible the process is. The quantity S expt (0) can be 
determined calorimetrically by performing a cooling ex- 
periment. We take Tq to be the melting temperature 
Tom, and uniquely determine the entropy of the super- 
cooled liquid at Tom by adding the entropy of melting 
to the crystal entropy Scr(Tim) at Tom- The latter is 
obtained in a unique manner by integration along a re- 
versible path from Tq — to Tq = Tim: 

Tom 

Scr(Tqm) = Scr(O) + J Cp t CB,dTo/To, 



here, S'cr(O) is the entropy of the crystal at absolute 
zero, which is traditionally taken to be zero in accordance 
with the third law, and C^crC^o) is the isobaric heat 
capacity of the crystal. This then uniquely determines 
the entropy of the liquid to be used in the right hand 
side in Eq. ((6]). We will assume that Scr(O) = 0. Thus, 
the experimental determination of 5 r cxp t(0) is required to 
give the lower bound to the residual entropy in Eq. ([2"]l, 
Experiment evidence for a non-zero value of S UmtfO) i s 
abundant as discussed by several authors |3.l5L llll4l^l23l : 
a textbook J38| also discusses this issue. Goldstein [n[ 
gives a value of Sr — 15.1 J/K mol for o-terphenyl from 
the value of its entropy at T = 2 K. We have given 
above a mathematical justification of S l eX pt(0) > in Eq. 
(|7|). The strict inequality proves immediately that the 
residual entropy cannot vanish for glasses, which justifies 
the curve Glass 1 in Fig. [2j 

The inequality in Eq. ([6]) takes into account any 
amount of irreversibility during vitrification; it is no 
longer limited to only small contributions of the order 
of 2% considered by several others [H (3 [3, IH, S3] , 
which makes our derivation very general. 

By considering the state Ao above to be a state Ao of 
the glass in a medium at some arbitrary temperature Tq 
below T 0g , we can get a generalization of Eq. ([6]): 



While we have only demonstrated the forward inequali- 
ties, the excess Sr— iS , exp t(0]can be computed in nonequi- 
librium thermodynamics [23, [30, |3fl 50, EH , which pro- 
vides a clear prescription for calculating the irreversible 
entropy generation. We do not do this here as we are 
only interested in general results, while the calculation of 
irreversible entropy generation will, of course, be system- 
dependent and will require detailed information. Gutzow 
and Schmelzer[12j| provide such a procedure with a sin- 
gle internal variable but under the assumption of equal 
temperature and pressure for the glass and the medium. 
However, while they comment that diS > whose eval- 
uation requires system-dependent properties, their main 
interest is to only show that it is negligible compared to 
d c S. 

We have proved Theorems [1] and [5] by considering only 
the system without paying any attention to the medium. 
For Theorem [TJ we require the second law, i.e. Eq. (QJ. 
This is also true of Eq. (0). The proof of Theorem [5] 
requires the constraint r b s < T rc i ax (To) for any To < Xbg, 
which leads to a nonequilibrium state. The same is also 
true of Eq. fllOj). 



s(r ) > SoxptCz^) 



S(T ) + JCpdTo/To. (8) 

To 



We again wish to remind the reader that all quantities 
depend on the path A— ^A , which we have not exhibited. 
By replacing To by the melting temperature Tom and T$ 
by To, and adding the entropy S(Tom) of the medium on 
both sides in the above inequality, and rearranging terms, 
we obtain (with S'l(Tom) = <Sscl(Tom) for the liquid) 

Sl(T m) + S(T 0M ) < S(T ) + S(T 0M ) - J C P dT Q /T , 

Tom 

(9) 
where we have also included the equality for a reversible 
process. This provides us with an independent derivation 
of the inequality given by Sethna and coworker [171 ] . 

It is also clear from the derivation of Eq. ([7]) that the 
inequality can be generalized to any temperature To < 
To g with the result 



<Sexpt(To) > Sscl(To), 



(10) 



with S , eX pt(T ) -> Sscl(To) as T -> T 0g from below. 
Thus, <S , oxp t(To) appears to have a form similar to that 
of Glassl in Fig. [3] but strictly lying below it. We have 
now verified the first statement in the Introduction. 



III. CONCLUSIONS 

We have considered the role of irreversible entropy gen- 
eration during isobaric vitrification to rigorously justify 
the two statements in the Introduction. They are valid 
regardless of how far the system is out of equilibrium. 
Thus, our results are very general and are not restricted 
by the small amount of irreversibility that is normally 
considered in the literature. The first statement shows 
that the instantaneous entropy S(To,t) must always be 
higher than S l exp t(To), which in turn must always be 
higher than Sscl(Io) of the equilibrated supercooled en- 
tropy. The second statement shows that the extrapola- 
tion of the calorimetrically measured entropy to absolute 
zero forms a strict lower bound to the residual entropy 
5r. As the former is usually positive, this proves that the 
residual entropy has to be at least as large as this value. 
From the first statement, it also follows that Glass2 is 
not realistic. 

The statements follow from considering the thermody- 
namic entropy that appears in the second law, and their 
validity is not affected by which equivalent statistical def- 
inition of entropy one may wish to use for the thermody- 
namic entropy, an issue that has been investigated by us 
recently [37j . 
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